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BENDING OF A RAY IN A RANDOM INHOMOGENEOUS MEDIUM! 


R.M. Jones 


Abstract 


When.a wave propagates through a random, inhomogeneous medium, 
the random component of the refractive index affects not only the 
random component of the ray direction, but also the mean component. 
Likewise, the mean component of the refractive index affects not 
only the mean component of the ray direction, but also the random 
component. This report derives equations based on a perturbation 
expansion that can be used to calculate the mean and random component 
of the ray direction in a medium in which the mean and random com- 
ponent of the refractive index vary arbitrarily in three dimensions. 
The equations are suitable for numerical integration in a ray tracing 
program. 


1. Summary 


This report develops a formula for the average ray direction and the mean 
square deviation from the average for propagation of a ray in a random, inhomo- 
geneous medium. Evaluation of the formula involves integrating three simultaneous 
equations along the average ray path and is thus suited for calculation by a ray 


tracing program. It is valid under the following restrictions: 
he Ray theory applies. 


a The medium is isotropic, so that the ray direction is the same 


as the wave normal direction. 


di The change in the average refractive index within a correlation 


length is negligible. 


4. The random component of the refractive index is much smaller than 


the average refractive index. 


„tú The mean sguare deviation of the ray direction from the average 


remains small. 
6. Three-point and higher correlations are neglected. 
Fa Moments of higher order than two are neglected. 


8. The ray direction is a Gaussian-distributed random variable. 


Presented at the U.R.S.I. International Symposium on Electromagnetic 
Wave Theory, USSR, Tbilisi, September 9-15, 1971. 


The formula actually gives the ensemble average < U> of the 
unit vector © in the ray direction. This ensemble average < o> 
determines the average ray direction. In addition < 0 > determines the 
mean square angular deviation < w“> from the average ray direction 


because <0 > =< coso > = exp( - >< o“ >). That is, 
< z = - 2 = oa 
d > Log < 0 > Log, ( < 0,7 «9,7 ) 


As in allthe equations in this paper, summation from 1 to 3 is implied 
over repeated indexes. 


The three components dg i of < 0 > are found by integrating the 


following three (for i= 1, 2, 3) differential equations along the average 


ray path: 
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where the total refractive index is 
n=< n> tu , 


the random part has a mean of zero with the following two -point cor- 


relation function 
«yu (Fa)> z<u? 2 C Un - Ta) , 


the magnitude of the average of the unit vector in the ray direction is 


sa» else > t= |< 0,><9,> ; 
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The derivation begins with a standard differential equation for 
the bending of a ray in an isotropic medium. As might be suspected 
from some of the limitations required for a valid result, the 
derivation involves expanding the ray equation in powers of i -< 9, > 


and keeping only linear terms. The resulting equation is solved 


43 


by successive approximations to give multiple integrals. Finally, an 
ensemble average is taken, neglecting three-point (and higher) correla- 


tions and moments of higher order than two. 


do 


2. Introduction 


In a homogeneous, random medium, the average ray direction remains 
constant, while the mean square deviation from the average increases with 
propagation distance through the medium. Formulas that give the mean 
square deviation of the ray direction from the average for small devia- 


tions are well known (e.g., Chernov, 1960). 


In an inhomogeneous medium the average ray direction is not constant 
and in fact depends on the randomness of the medium. The purpose of this 
report is to calculate the average ray direction and the mean square 
deviation from the average for a ray traveling through an inhomogeneous 
medium. The calculations for an inhomogeneous medium are much more com- 
plicated than for a homogeneous medium. They would be even more com- 
plicated without limiting the applicability of the calculations to situa- 
tions that satisfy the following restrictions. 


M Ray theory applies. 


ka The medium is isotropic, so that the ray direction is the 


same as the wave-normal direction. 


Je The change in the average refractive index within a correlation 
length is negligible. This is necessary because the derivation depends 
on dividing the medium into segments small enough that the gradient of 
the average refractive index is nearly constant but still much larger 
than the correlation length for the random component of the refractive 


index. 


CMS The random component of the refractive index is much smaller 
than the average refractive index. This is necessary because the deri- 
vation involves an expansion in powers of the ratio of the random com- 
ponent of the refractive index to the average refractive index. Most 


media probably satisfy this restriction. 


ins s 


5. lhe mean square deviation of the ray direction from the 
average ray direction remains small. This is necessary because the 
derivation involves an expansion in powers of the ratio of difference 
between the instantaneous unit vector in the ray direction and the aver- 
age of a unit vector in the ray direction to the average of a unit vector 
in the ray direction. This approximation will eventually break down, 
even in a homogeneous medium, for long path lengths. The final solu- 
tion is therefore valid only for relatively short path lengths. 

6. Three-point and higher correlations are neglected. This 
limitation simplifies the calculations. Higher correlations could be 
added although their effect would probably be small in most cases. 

T Moments of higher order than two are neglected. This 
simplifies the calculations. Higher moments could be added. 

8. The ray direction is a Gaussian-distributed random 


variable. 


3. Derivation Method 


The average ray direction and the mean square deviation from 
the average are determined by the average of a unit vector in the ray 
direction. It is clear that the average of a unit vector in the ray direc- 
tion gives the average ray direction. It is not immediately clear that 
it also gives the mean square deviation from the average. Section 4 
shows this to be true, if the ray direction is a Gaussian-distributed 
random variable, and then derives the required formula. | 

Lhe calculation of the average of a unit vector in the ray direc- 
tion for an inhomogeneous medium can be broken into two parts. First, 
we consider a medium in which the average refractive index has a con- 
stant gradient and the correlation functions for the random part of the 


refractive index are constant. Second, we simply consider the 


T - 


inhomogeneous medium to be made up of segments that satisfy the above 
criteria, and we then write the final solution as an integral over these 
segments. 

We begin by writing the general ray equation for the bending of 
a ray in an isotropic medium in terms of a unit vector o in the ray 
direction. We then choose a coordinate system in which the gradient of 
the average refractive index is in the z direction, and we rewrite the ray 
equation so that z is the independent variable instead of ray path length. 
We let our medium be a slab bounded by z= 0 and z 7 Az. 

| Next we convert to tensor notation for ease in factoring vectors 
out of dot products. We then write the unit vector. © às the sum of the 
average unit vector < o> plus a small perturbation AG and expand the 
ray equation in powers of Nc . We then neglect all powers of Ko but 
„linear terms. Keeping higher powers would give terms in the final re- 
sult containing three-point and higher correlations or terms proportional 
to the square of the slab thickness Az. 

We now write the refractive index as.a sum of its average value 
< n> plus its random component NS. The average refractive index 
varies linearly through the slab because we assume a constant gradient 
of <n> inthe slab. Keeping the linear variation of < n> in the 
equation gives terms in the final result proportional to the square of 
the slab thickness. Since we want ultimately to neglect such terms, we 
can neglect them now by assuming < n> constant within the slab. We 
then expand the equation in a power series of u/ < n> and discard 
powers larger than 2. This corresponds to neglecting moments of 
higher order than 2 in the ensemble average. 

We now solve the differential equation by successive approxima- 
tions, beginning with No 20. This gives multiple integrals over the 


slab thickness. We neglect triple integrals and higher order integrals 


because they lead to terms in the final result that contain three-point 
and higher order correlations or terms proportional to the square and 
higher powers of the slab thickness. We also discard other terms in 
the successive approximation solution proportional to the square of the 
slab thickness. Taking the ensemble average of the remaining single 
and double integrals gives the change in the average of the unit vector 
after having gone through the slab. 

lhe change in the average value of the unit vector is proportional 
to the approximate straight-line distance traveled in the slab. Taking 
the limit of small slab thickness allows rewriting the equation as a 
differential equation for the average of the unit vector in the ray direc- 
tion. We assume that this differential equation applies at each point in 
the inhomogeneous medium. 

Finally, the special case of a homogeneous medium is compared 


with the result of Chernov (1960). 


4, The Derivation 


We assume that the refractive index n is composed of an aver- 


age value plus a random component 
n=<n> tu (1) 


where the brackets < and > indicate an ensemble average. The ray 


equation in an isotropic medium can be written 


dog) - 
A - Vn (2) 


e 
where O is a unit vector in the ray direction (same as the wave normal 
direction since the medium is isotropic), and s is geometrical path 

EX 
length. We are considering a medium in which V < n» and all the 


statistical properties of u are constant. We choose the z axis parallel 


mo 


with V « n». Our medium is a slab with thickness Az and with 
boundaries perpendicular to the z axis. We assume that a representa- 
tive ray enters the slabat z - 0 and leaves at z=Az. Thus no 
reflection occurs within the slab, only a small amount of bending. We 
now change the independent variable in the ray equation from s to z 


by means of the following equation: 


ds lv<n>| 
E c (3) 
^X MR us EE > 
Using 
Un "m ! 
da 9 * (4) 


along with (3) in (2) gives 


n |V € n>] e eh ae — 
Se (vn -5 (G+ Fn) ) i (5) 


Equation (5) can be rewritten in tensor notation 


do: stre a>l V.n-0,(0,7.n) lv.¢n >| 6 479,9, Vn 
ER oH S e n 


oe — 
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E (6) 
dz n g. v n> n Z v, n 


where summation from ] to 3 is implied over repeated indices, and 
al is the Kronecker ô . Equation (6) explicitly gives each component 


of (5) for i equal 1, 2, ae X. 


u | 
O. o> tho, i (7) 


Then 
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Using (8) and (9) in (6) gives 


do, lv<n>| V.H fa -< Um c9.» ww G7 h0.-«G > bo, - DO. AO, ) 
A od J 1J 1 J J 1 1 J 1 J 
WF eae £ 
dz c eN kn . (14 Vk n> ay 
k 


<o> | 
(10) 


We can now expand (10) in powers of Ao, discard powers of 2 and higher, 


and collect terms, changing dummy indices where necessary to give 


do, lv<n>| V... JA 
—7 ————— —— (5..-<a.><0.> + 
"M ites VA > “J 7 J 


+( -<0,>0, "cd b. (5..-<0.><0.>) 
l- X 1 vs M ij i j p 


Substituting (1) into (11), we obtain 


do, b<o>rliv<ar| V; en»tV.u 
— —RBHO i (and, he, ) (12) 
az <0> V<n> 1+ + S MN UN 

«n» 


2105 


where 


l 
À.. 6 «0 eG. > (13) 
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Expanding (12) in powers of ų and discarding powers higher than 2 gives 


do; I«o»||v« n7] 


e s 
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(15) 
which can be rewritten 
do, 
Te O. + Ba. AO (16) 


where 
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(18) 
Equation (16) can be rewritten in integral form 
Z Z 
AG (z) = | a (zı) dz, + | 84,0) AG, (zi) dz4 : (19) 
O O 


a1. 


Equation (19) can be solved by successive approximations. For the 


zeroth order approximation, we choose 


^ 0. (z =0 . (20) 


Substituting (20) into (19) gives for the first-order approximation 


Z 
A 0. (2) = | U. (zı) dz, > (21) 


O 


Substituting (21) into (19) gives the second-order approximation 


Az 
Ao, (A z) = | a. (zı) dz, + 
o (22) 
Az Z1 
* | 8i, (21) | k (Za) dza dz, 
o o) 


It is unnecessary to go to higher order approximations, because that 
would add only terms with triple and higher order integrals and terms 
provortional to the sguare of the slab ťhickness. 

Now we take the ensemble average of (22). If we substitute (17) 
into the single integral in (22) and expand the product, we have 6 terms 
in the single integral. Appendix A evaluates the ensemble averages of 
these single integrals. If we substitute (17) and (18) into the double 
integral in (22) and expand the products, we have 35 terms in the double 
integral. Some of these terms can be separated into the product of two 
single integrals and evaluated using Appendix A. Many of the ensemble 
averages of the remaining double integrals are zero because of the 
assumptions we made. Appendix B evaluates the ensemble averages of 


the non-zero double integrals. 


lx. 


Substituting the formulas for the integrals from Appendices A 
and B into (22) and discarding terms proportional to the square of the 


slab thickness gives 


< A0 > =A, Y< n> ín T ) 58+ 


(23) 
< < - < < ^ 
TB. Sk ^ kt Cy, n>V, n> "rd n» t FE n» Ki, )A S 
The integrals I, J; » and Kij are given by (B-16), (B-20), and 
(B-21). Notice that we have changed from slab thickness to average 
path length in the slab by using (3). We can now let the slab thickness 
approach zero, and rewrite (23) as a differential equation fọr the aver- 


age of the unit vector in the ray direction. 


dI». 


i 
ise au 9 2 ( £t Suz) 
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Substituting (13) and (14) into (24) gives (after much algebra) 
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where 
aoe ey <0 >| = | <9,><0,> n (26) 
da EA nzie< APT. (27) 


the integral I is given by (B-16), Ji is given by (B -20), Kt is given 
by (B-21), and om that appears in (B-16), (B-20), and (B-21) is the 


two-point correlation function defined by (B-4). 


5. Significance of the Result 


Equation (25) gives a system of 3 (for i= 1, 2, 3) differential 
equations that can be numerically integrated along the average ray 
path to give the final value of < 0> for a given set of initial conditions 
at the beginning of the ray path. Equation (25) is in a form that can be 
integrated by a ray tracing program. 

The average of the unit vector in the ray direction < o> deter- 
mines the average ray direction. In addition, < c > determines the 
mean square angular deviation < w“ > from the average ray direction 


in the following way: 


on 


Let © be the angle between the instantaneous ray direction and 


the average ray direction. Then 


0. 50,» 
doti. vot TC (28) 
Taking the ensemble average of (28) gives 
úd ZAR x 2 , 
< 0 > 
«cog > = pe = cas 3« 9g: e (29) 


The angle © is a random variable with mean zero. If in addition it is 


Gaussian distributed, then 
<cos > = exp( -< ep“ >/2) i (30) 
| Combining (29) iná (30) gives 
KY“ > =.-2 log: 0 ok chop e Go s (31) 
6. Comparison with Previous Results for a Homogeneous Medium 
For the special case that 


T M L0. (32) 


(25) shows that « o> does not change in direction but only in magnitude, 


and we have 


4502. db E 4 sále ON. Aou. EU 
S re) J 


If in addition Kg is constant, then (33) can be integrated to give 


Boe 
<n>? 


«g»? 2«g >? +2 S (34) 
S O 
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Combining (31) and (34) gives for the mean square angular deviation 
from the average ray direction after traveling a distance s through 


the medium 


Bee 


< o“ T2 - log, (expl -< v? > )*2 > noe s) (35) 


where < v“ 7, is the initial mean square angular deviation from the 
average ray direction. For the special case of a Gaussian correlation 


function 


C (r) = exp (-r? / a^) (36) 


(35) becomes 
« © <a" log, ( exp( « © T 4 2 JI] >; (37) 


If we now take the case for 


2 n. 
< © aj = 0 (38) 
and 
<a> s] (39) 
then (37) becomes 
2 p di 2 5 
< © > a = log. (1 acu" Jm D i (40) 


For s small enough, we can take the first term in the power series 


for (40) to give 


<n? > m 4«y?7 A= | (41) 


which agrees with the result given by Chernov (1960, page 17). 
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APPENDIX A 


Evaluation of the Ensemble Averages of the Single Integrals 


Substituting (17) into the single integral in (22) gives 6 integrals 


whose ensemble averages must be calculated. 
The first of these integrals is not a random variable and can 


easily be evaluated and then rewritten using (3). 


Az <o> . Yxn 
(A-1) 


de = dg arnesa 66 
| «o»||v <n>| 


The ensemble average of the second integral is zero, since y 


has an average of zero. 
Az TL ^z 
| uaz -|«u»dz-| 0dz-0 {A -2) 
O O O 


The ensemble average of the third integral is easily calculated 


and then rewritten using (3) 
Az <0> + J< n> 


€u?» da =<? > hy @ <*> ^s 


I«o»||v«n7| 


(A-3) 


The ensemble average of the fourth integral is zero because y 


has an average of zero. 


Az Az , Az 
| Vu dz | ¥<p> dz z (A -4) 
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The ensemble average of the fifth integral is zero because we 


have assumed that < u? > is constant within the slab. 


aa 


Az. i GEM. | Az 
[uF uaz = > | <¥ (u?)> dz = | ¥<p?> dz = (A - 5) 


N |1 


2 E 
O O 


The ensemble average of the sixth integral is zero because we 


neglect third moments. 


^z "P : ^z 
|u? Vu dz -2|V«u*»dz-0 (A - 6) 
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APPENDIX B 


Evaluation of the Ensemble Averages of the Double Integrals 


Because we are neglecting three-point and higher correlations, 
only four of the 36 double integrals do not factor into the product of two 
single integrals or have ensemble averages that are not zero. 


The first of these integrals is 


hz Zi 
i u (zı) I" (za) dz. dz, , (B -1) 


O O 


which can be rewritten 


Az p 21 
< jH (zı) u (z2) dz, D : (B -2) 


O O 
ľaking the ensemble average inside the integral gives 


Am 74 Az Z4 
| «cu e (z. - Zi) dz. dz, x <ti“ > | | C (z2 - Z1) dzez dz: (B -3) 


O O O O 


where the correlation function Ca is defined by 
<u (ridu(rg)> "59^ 2 C (m, - ra) . (B -4) 


We now make a change of integration variables by letting 


"E IL (B -5) 
and 
Z 7 Zg - Zi ; (B -6) 


„20 = 


Then we have 


dza dz, = dz dz, , (B -7) 
and (B-3) becomes 
^z 
ae ^z o 
«xu? | | C (z) dz dz. +<y*> | | C (z) de dm. .. (B -8) 
: x T i a: | : 
^ 20 — 2(zg - Nz) 


We assume that the correlation length is much less than the slab thick- 
ness. Thus, the correlation function E is nearly zero except for z 
near zero. Inthat case, we can extend the lower limit of integration 
of the inner integrals in (B-8) to minus infinity without adding a signifi- 


cant amount to the integrals. Therefore, (B-8) becomes 
^z , 
2 o Az 9 
<u “> | | C (z) dz dz, +< u^» | | C (z) dz dza . (B -9) 


O - 00 ^z - 00 
č 
Combining the integrals in (B-9) gives 


As o 
zy? | | C (z) dz dz. . (B -10) 


O = CO 


Since the inner integral is independent of z,, the outer integral can 


easily be evaluated so that (B-10) becomes 


O 


F | C (2) dz . (B-11) 


an CO 
We assume that the correlation function satisfies 


es (-z) = C (z) : (B-12) 


„si. 


Using (B-12), (B-11) becomes 
KU? » Ae | c (z) A | (B-13) 


O 


We want to express (B-13) in terms of the average path length of the ray 
in the slab rather than the slab thickness ^z. We can use (3) to make 
this change. In addition, it is useful to change the integration variable 
in (B-13) from distance into the slab to distance traveled along the aver- 
age ray path. Equation (3) can also be used to make this change. Also, 
the correlation function M can be anisotropic so that it depends on the 
ray direction. This dependence should now be shown explicitly. Thus, 


(B-13) becomes 


à eo) ^ oily 2 E ved 
uj > ——-—- 0 às | C (op = ) ap „A BUT 
|< o>||v< n>| - s codi 


lhe final result, giving the equivalence of (B-1) and (B-14), is 


Nz Za ; eg sa 2 
«gue 

| ute.) | u (z2) dza diii eB dun oa I^s (B21) 

|<o>||v<n>| 
O O 
where 

c RE > S p <I> 
E cor € dao. eee. 


Similarly, the ensemble averages of the other three non-zero integrals 


can be calculated to give 


me be a 


^z Z4 


2 
| u (zi) | Vy Hz) dz. dz, =< n> C P ise J As 
: f |<o>||v<n>| 


(B-17) 
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|<o>||v<n>| 


O O 
(B -18) 
and 
^z Z1 
ANTO, [v Hag) dz, dz; ) = (sm zar S ET kt? 
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(B -19) 
where 
iR ^T» E " Ca > ) i 
J —— V. Oi (P Ve do (B -20) 
O 
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[e] 
SN 2 
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